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(54) Title: LAGUERRE FONCTION FOR AUDIO CODING 




(57) Abstract: Coding (1) of an 
audio signal is provided including 
estimating (110) a position of a 
transient signal component in the 
audio signal, matching (111, 112) a 
shape function on the transient signal 
component in case the transient signal 
component is gradually declining 
after an initial increase, which shape 
function has a substantially exponential 
initial behavior and a substantially 
logarithmic declining behavior; and 
including (15) the position and shape 
parameters describing the shape 
function in an audio stream (AS). 



WO 01/69593 PCT/EP01/02424 

1 



LAGUERRE FUNCTION FOR AUDIO CODING 

The invention relates to coding of audio signals, in which transient signal 
components are coded. 

The invention further relates to decoding of audio signals. 
The invention also relates to an audio coder, an audio player, an audio system, 
5 an audio stream and a storage medium. 

The article from Purnhagen and Edler, "Objektbasierter Analyse/Synthese 
Audio Coder fur sehr niedrige Datenraten", ITG Fachbericht 1998, No. 146, pp. 35-40 
discloses a device for coding of audio signals at low bit-rates. A model-based Analysis- 

10 Synthesis arrangement is used, in which an input signal is divided in three parts: single 

sinusoids, harmonic tones, and noise. The input signal is further divided in fixed frames of 32 
ms. For all blocks and signal parts, parameters are derived based on a source-model. To 
improve the representation of transient signal parts, an envelope function a{t) is derived from 
the input signal and applied on selected sinusoids.. The envelope function consists of two line 

1 5 segments determined by the parameters r at k> /'dec, 'max as shown in Fig. 1 . 

An object of the invention is to provide audio coding that is advantageous in 
terms of bit-rate and perception. To this end, the invention provides a method of coding and 
decoding, an audio coder, an audio player, an audio system, an audio stream and a storage 
20 medium as defined in the independent claims. Advantageous embodiments are defined in the 
dependent claims. 

A first embodiment of the invention comprises estimating a position of a 
transient signal component in the audio signal, matching a shape function on the transient 
signal component in case the transient signal component is gradually declining after an initial 
25 increase, which shape function has a substantially exponential initial behavior and a 

substantially logarithmic declining behavior; and including the position and parameters 
describing the shape function in an audio stream. Such a function has an initial behavior 
substantially according to t" and a declining behavior after the initial increase substantially 
according to where /is a time, and n and a are parameters which describe a form of the 
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shape function. The invention is based on the insight that such a function gives a better 
representation of transient signal components while the function may be described by a small 
number of parameters, which is advantageous in terms of bit-rate and perceptual quality. The 
invention is especially advantageous in embodiments where transient signal components are 
5 separately encoded from a sustained signal component, because especially in these 
embodiments a good representation of the transient signal components is important. 

According to a further aspect of the invention, the shape function is a Laguerre 
function, which is in continuous time given by 
c-re-" (1) 

10 where c is a scaling parameter (which may be taken one). In a practical embodiment, a time- 
discrete Laguerre function is used. 

Transient signal components are conceivable as a sudden change in power (or 
amplitude) level or as a sudden change in waveform pattern. Detection of transient signal 
components as such, is known in the art. For example, in J. Kliewer and A. Mertins, ' Audio 

1 5 subband coding with improved representation of transient signal segments', Proc. of 

EUSIPCO-98, Signal Processing IX, Theories and applications, Rhodos, Greece, Sept. 1998, 
pp. 2345-2348, a transient detection mechanism is proposed, that is based on the difference in 
energy levels before and after an attack start position. In a practical embodiment according to 
the invention, sudden changes in amplitude level are considered. 

20 In a preferred embodiment of the invention, the shape function is a generalized 

discrete Laguerre function. Meixner and Meixner-like functions are practical in use and give 
a surprisingly good result. Such functions are discussed in A.C. den Brinker, 'Meixner-like 
functions having a rational z-transform', Int. J. Circuit Theory AppL, 23, 1995, pp. 237-246. 
Parameters of these shape functions are derived in a simple way. 

25 In another embodiment of the invention, the shape parameters include a step 

indication in case the transient signal component is a step-like change in amplitude. The 
signal after the step-like change is advantageously coded in sustained coders. 

In another preferred embodiment of the invention, the position of the transient 
signal component is a start position. It is convenient to give the start position of the transient 

30 signal component for adaptive framing, wherein a frame starts at the start position of a 
transient signal component. The start position is used both for the shape function and the 
adaptive framing, which results in efficient coding. If the start position is given, it is not 
necessary to determine the start position by combining two parameters as would be necessary 
in the embodiment described by Edler. 
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The aforementioned and other aspects of the invention will be apparent from 
and elucidated with reference to the embodiments described hereinafter. 



In the drawings: 

5 Fig. 1 shows a known envelope function, as already discussed; 

Fig. 2 shows an embodiment of an audio coder according to the invention; 
Fig. 3 shows an example of a shape function according to the invention; 
Fig. 4 shows a diagram of first and second order running central moments of 
an input audio signal; 

10 Fig. 5 shows an example of a shape function derived for an input audio signal; 

Fig. 6 shows an embodiment of an audio player according to the invention; 

and 

Fig. 7 shows a system comprising an audio coder and an audio player; 

The drawings only show those elements that are necessary to understand the 

1 5 invention. 



Fig. 2 shows an audio coder 1 according to the invention, comprising an input 
unit 10 for obtaining an input audio signal x(Y). The audio coder 1 separates the input signal 
into three components: transient signal components, sustained deterministic components, and 

20 sustained stochastic components. The audio coder 1 comprises a transient coder 1 1, a 
sinusoidal coder 13 and a noise coder 14. The audio coder optionally comprises a gain 
compression mechanism (GC) 12. 

In this advantageous embodiment of the invention, transient coding is 
performed before sustained coding. This is advantageous because transient signal 

25 components are not efficiently and optimally coded in sustained coders. If sustained coders 
are used to code transient signal components, a lot of coding effort is necessary, e.g. one can 
imagine that it is difficult to code a transient signal component with only sustained sinusoids. 
Therefore, the removal of transient signal components from the audio signal to be coded 
before sustained coding is advantageous. A transient start position derived in the transient 

30 coder is used in the sustained coders for adaptive segmentation (adaptive framing) which 
results in a further improvement of performance of the sustained coding. 

The transient coder 1 1 comprises a transient detector (TD) 1 10, a transient 
analyzer (TA) 1 1 1 and a transient synthesizer (TS) 112. First, the signal x(t) enters the 
transient detector 110. This detector 1 10 estimates if there is a transient signal component, 
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and at which position. This information is fed to the transient analyzer 111. This information 
may also be used in the sinusoidal coder 13 and the noise coder 14 to obtain advantageous 
signal-induced segmentation. If the position of the transient signal component is determined, 
the transient analyzer 111 tries to extract (the main part of) the transient signal component. It 
5 matches a shape function to a signal segment preferably starting at an estimated start 

position, and determines content underneath the shape function, e.g. a (small) number of 
sinusoidal components. This information is contained in the transient code Cj. The transient 
code Ct is furnished to the transient synthesizer 112. The synthesized transient signal 
component is subtracted from the input signal x(t) in subtractor 16, resulting in a signal x\. In 

10 case, the GC 12 is omitted, x\ — X2- The signal xi is furnished to the sinusoidal coder 1 3 where 
it is analyzed in a sinusoidal analyzer (SA) 130, which determines the (deterministic) 
sinusoidal components. This information is contained in the sinusoidal code Cs. From the 
sinusoidal code Cs, the sinusoidal signal component is reconstructed by a sinusoidal 
synthesizer (SS) 131. This signal is subtracted in subtractor 17 from the input x 2 to the 

15 sinusoidal coder 13, resulting in a remaining signal X3 devoid of (large) transient signal 

components and (main) deterministic sinusoidal components. Therefore, the remaining signal 
X3 is assumed to mainly consist of noise. It is analyzed for its power content according to an 
ERB scale in a noise analyzer (NA) 14. The noise analyzer 14 produces a noise code Cn. 
Similar to the situation in the sinusoidal coder 13, the noise analyzer 14 may also use the start 

20 position of the transients signal component as a position for starting a new analysis block. 
The segment sizes of the sinusoidal analyzer 130 and the noise analyzer 14 are not 
necessarily equal. In a multiplexer 15, an audio stream AS is constituted which includes the 
codes Cj, Cs and Cn- The audio stream AS is furnished to e.g. a data bus, an antenna system, 
a storage medium etc. 

25 In the following, a representation of transient signal components according to 

the invention will be discussed. In this embodiment, the code for transient components Cj 
consists of either a parametric shape plus the additional main frequency components (or other 
content) underneath the shape or a code for identifying a step-like change. According to a 
preferred embodiment of the invention, the shape function for a transient that is gradually 

30 declining after an initial increase, is preferably a generalized discrete Laguerre function. For 
other types of transient signal components, other functions may be used. 

An example of a generalized discrete Laguerre function, is a Meixner function. 
A discrete zeroth-order Meixner function g{t) is given by: 
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g(0 = ^(W 2 ) i/2 £' (2) 

where / = 0,1,2,.. and (b) t = 6(6+1).. .(6+M) is a Pochhammer symbol. The parameter b 
denotes an order of generalization (b > 0) and determines the initial shape of the function: 
approximately / oc t (b ~ l)/1 for small t. The parameter £ denotes a pole with 0 < 1 and 
5 determines the decay for larger t. The function g(f) is a positive function for all values of t. 
For b = 1, a discrete Laguerre function is obtained. Furthermore, for b = 1, the z-transform of 
g is a rational function in z and can thus be realized as an impulse response of a first order 
infinite impulse response (IIR) filter. For all other values of b there is no rational z- 

OO 

transform. The function g(t) is energy normalized, i.e. ^ g 2 (/) = 1 . The zeroth-order 

10 Meixner-function may be created recursively by: 
g(0) = (l-Z 2 ) b/2 (3) 

g(l) = ^~^^g{t - 1) for t >0 (4) 

In another embodiment according to the invention, Meixner-like functions are 
used, because they have a rational z-transform. An example of a Meixner-like function is 
15 shown in Fig. 3. A discrete zeroth-order Meixner-like function h(i) is given by its z- 
transform: 



H(z) = C e 



Z 



where a = 0,1,2,... and C a is given by: 



(5) 



✓1 e2xa+l/2 _ er2x(a+l)/2 

C a = \ I J . - = \ I j = (6) 



n=0 



^ 2 



^ 2 N 



20 where P a is an ath order Legendre polynomial, given by: 

The parameter a denotes the order of generalization (a is a non-negative integer) and £ is the 
pole with 0 < £ < 1 . The parameter a determines the initial shape of the function: / cct a for 
small (. The parameter £ determines the decay for large t. The function h is a positive 
25 function for all values of t and is energy normalized. For all values of a, the function h has a 
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rational z-transform and can be realized as the impulse response of an IIR filter (of order 
a+1). 

The function h(t) can be expressed in a finite discrete Laguerre-series 

according to: 

> xo=2X<uo (8) 

where ^ are discrete Laguerre functions, see the article of A.C. den Brinker. B m is given by: 

f „\ 

(9) 



B =C ^ 



2\a+l/2 



a 



(i-rr 

First and second order running central moments of a given function f(t) are 



defined by: 



Z: e-w 2 (» 

io W)= — do) 



E- C-*o-r,(*)) 2 / 2 (0 
r,(Jt)= r (ii) 

11 z:::,/ 2 (o 

where ko is the start position of the transient signal component. 

With a good estimation of the running moments T\ and T2 of an input audio 
signal (take f(t) = x(t) in equations 10 and 1 1), the shape parameters may be deduced. 

1 5 Unfortunately, in real data a transient signal component is usually followed by a sustained 
excitation phase, disturbing a possible measurement of the running moments. Fig. 4 shows 
the first and second order running central moments of an input audio signal. It appears that 
the running moments initially increase linearly from the assumed starting position and later 
on tend to saturate. Although the shape parameters may be deduced from this curve, because 

20 the saturation is not as clear as desired for parameter extraction, i.e. it is not clear enough at 
which k good estimates of Tj and T 2 are obtained. In an advantageous embodiment of the 
invention, a ratio in initial increase of the running moments T } and T2 is used to deduct the 
shape parameters. This measurement is advantageous in determining b (and in case of the 
zeroth-order Meixner function a), since b determines the initial behavior of the shape. From a 

25 ratio between slopes of running moments T\ and T2 a good estimation for b is obtained. From 
simulation results has been obtained that to a very good degree, a linear relation exists 
between the ratio slope Tj/ slope T 2 and the parameter b, which is, in contrast to a Laguerre 
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function, slightly dependent on the decay parameter £ As a description may be used (derived 
by experiments): 

for Meixner: slope Tj/ slope T 2 = b + 1/2 (12) 

for Meixner-like: s7o/?e Tj/slope T 2 = 2a + 5/2 (13) 

5 wherein a £ dependence is ignored. Because Tj and T 2 are zero for k = s/ope Tj/ slope T 2 
may be approximated by Tj/ T 2 for a suitable £. 

The pole £ of the shape may be estimated in the following way. A second 
order polynomial is fitted to a running central moment, e.g. Tj. This polynomial is fitted to a 
signal segment of T\ with observation time Tsuch that leveling off is clearly visible, i.e. a 
1 0 clear second order term in the polynomial fit at T. Next, the second-order polynomial is 

extrapolated to its maximum and this value is assumed to be the saturation level of T\. From 
this value for T\ and b, £is calculated with use of equations 2 and 10, with f(t) = g(t). For a 
Meixner-like function, £is calculated from the value for T\ and a, with use of equations 8-10, 
with/(?; = h(t). 

15 A procedure for estimation of the decay parameter g is as follows: 
start with some value of T 

fit a second order polynomial to the data on 0 to 7, i.e. T } (/) « c 0 +c, t + c 2 t 2 for t = [0,7] 
where co,/,2 are fitting parameters 

check if the quadratic term of this polynomial is essential at t = T: 
20 T x (T) < (1 - £)(c 0 + C \T) where 8 represents a relative contribution of the quadratic term at / = 

T. 

if this is satisfied, then extrapolate T\(t) to its maximum and equate this with T\\ 

2 



calculate the decay parameter % from T\ and b (or a) 
25 For Meixner-like functions, the shape parameter a is preferably rounded to integer values. 

Fig. 5 shows an example of a shape function derived for an input audio signal. 
Some pre-processing, like performing a Hilbert transform of the data, may be 
performed in order to get a first approximation of the shape, although pre-processing is not 
essential to the invention. 
30 When the value at which the running moments saturate is large, i.e. in the 

order of segment/ frame length, the Meixner (-like) shape is discarded. In case the transient is 
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a step-like change in amplitude, the position of the transient is retained for a proper 
segmentation in the sinusoidal coder and the noise code. 

After the start position and the shape of a transient have been determined, the 
signal content underneath the shape is estimated. A (small) number of sinusoids is estimated 
5 underneath the shape. This is done in an analysis-by-synthesis procedure as known in the art. 
The data that is used to estimate the sinusoids, is a segment which is windowed in order to 
encompass the transient but not any consequent sustained response. Therefore, a time 
window is applied to the data before entering the analysis-by-synthesis method. In essence, 
the signal which is considered extends from the start position to some sample where the 

10 shape is reduced to a certain percentage of its maximum. The windowed data may be 

transformed to a frequency domain, e.g. by a Discrete Fourier Transform (DFT). In order to 
avoid low-frequency components, which presumably extend beyond the estimated transient, a 
window in the frequency domain is also applied. Next the maximum response is determined 
and the frequency associated with this maximum response. The estimated shape is modulated 

15 by this frequency, and the best possible fit is made to the data according to some 

predetermined criterion, e.g. a psycho-acoustic model or in a least-squares sense. This 
estimated transient segment is subtracted from the original transient and the procedure is 
repeated until a maximum number of sinusoidal components is exceeded, or hardly any 
energy is left in the segment. In essence, a transient is represented by a sum of modulated 

20 Meixner functions. In a practical embodiment, 6 sinusoids are estimated. If the underlying 
content mainly contains noise, a noise estimation is used or arbitrary values are given for the 
frequencies of the sinusoids. 

The transient code Cj includes a start position of a transient and a type of 
transient. The code for a transient in the case of a Meixner (-like) shape includes: 

25 the start position of the transient 

an indication that the shape is a Meixner (-like) function 
shape parameters b (or a) and £ 

modulation terms: TVp frequency parameters and amplitudes for (co)sine modulated shape 
In case that the transient is essentially a sudden increase in amplitude level 
30 where there is no clear decay in this level (relatively) shortly after the starting position, the 
transient cannot be encoded with a Meixner (-like) shape. In that case, the start position is 
retained in order to obtain proper signal segmentation. The code for step-transients includes: 
the start position of the transient 
an indicator for the step 
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The performance of the subsequent sustained coding stages (sinusoidal and 
noise) is improved by using the transient position in the segmentation of the signal. The 
sinusoidal coder and the noise coder start at a new frame at the position of a detected 
transient. In this way, one prevents averaging over signal parts, which are known to exhibit 
5 non-stationary behavior. This implies that a segment in front of a transient segment has to be 
shortened, shifted or to be concatenated with a previous frame. 

The audio coder 1 according to the invention optionally comprises a gain- 
control element 12 in front of the sustained coders 13 and 14. It is advantageous for the 
sustained coders, to prevent changes in amplitude level. For a step-transient, this problem is 
1 0 solved by using a segmentation in accordance with the transients. For transients represented 
with an shape, the problem is partly solved by extracting the transient from the input signal. 
The remnant signal still may include a significant dynamic change in amplitude level, 
presumably shaped similar to the estimated shape. In order to flatten the remnant signal, the 
gain control element may be used. A compression rate may be defined by: 

15 gc(t) = (12) 

S W 1 + dhit) 

wherein h{f) is the estimated shape and d is a parameter describing a compression rate. The 
gain-control element assumes that after a transient, a stationary phase occurs with amplitude 
excursions amounting to about 0.2 times the maximum in the estimated shape. A ratio r is 
defined by: 

20 r = M r -0.2M e 

0.2M e 

wherein M T is the maximum of the remnant signal. 

The compression rate parameter d is equal to r if r > 2, otherwise d is taken 0. For the 
compression, only d needs to be transmitted. 

Fig. 6 shows an audio player 3 according to the invention. An audio stream 

25 AS', e.g. generated by an encoder according to Fig. 2, is obtained from a data bus, an antenna 
system, a storage medium etc. The audio stream AS is de-multiplexed in a de-multiplexer 30 
to obtain the codes C T \ Cs' and C N '. These codes are furnished to a transient synthesizer 31, 
a sinusoidal synthesizer 32 and a noise synthesizer 33 respectively. From the transient code 
Ct\ the transient signal components are calculated in the transient synthesizer 31. In case the 

30 transient code indicates an shape function, the shape is calculated based on the received 

parameters. Further, the shape content is calculated based on the frequencies and amplitudes 
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of the sinusoidal components. If the transient code Or' indicates a step, then no transient is 
calculated. The total transient signal yj is a sum of all transients. 

In case the decompression parameter d is used, i.e. if derived in the coder 1 
and included in the audio stream AS\ a decompression mechanism 34 is used. The gain 
5 signal g{t) is initialized at unity, and the total amplitude decompression factor is calculated as 
the product of all the different decompression factors. In case the transient is a step, no 
amplitude decompression factor is calculated. 

From two subsequent transient positions, a segmentation for the sinusoidal 
synthesis SS 32 and the noise synthesis NS 33 is calculated. The sinusoidal code C s is used to 
1 0 generate signal ys, described as a sum of sinusoids on a given segment. The noise code Cn is 
used to generate a noise signal ^n. Subsequent segments are added by, e.g. an overlap-add 
method. 

The total signal y(t) consists of the sum of the transient signal yj and the 
product of the amplitude decompression g and the sum of the sinusoidal signal ys and the 

1 5 noise signal y^. The audio player comprises two adders 36 and 37 to sum respective signals. 
The total signal is furnished to an output unit 35, which is e.g. a speaker. 

Fig. 7 shows an audio system according to the invention comprising an audio 
coder 1 as shown in Fig. 2 and an audio player 3 as shown in Fig. 6. Such a system offers 
playing and recording features. The audio stream AS is furnished from the audio coder to the 

20 audio player over a communication channel 2, which may be a wireless connection, a data 
bus or a storage medium. In case the communication channel 2 is a storage medium, the 
storage medium may be fixed in the system or may also be a removable disc, memory stick 
etc. The communication channel 2 may be part of the audio system, but will however often 
be outside the audio system. 

25 It should be noted that the above-mentioned embodiments illustrate rather than 

limit the invention, and that those skilled in the art will be able to design many alternative 
embodiments without departing from the scope of the appended claims. In the claims, any 
reference signs placed between parentheses shall not be construed as limiting the claim. The 
word 'comprising' does not exclude the presence of other elements or steps than those listed 

30 in a claim. The invention can be implemented by means of hardware comprising several 
distinct elements, and by means of a suitably programmed computer. In a device claim 
enumerating several means, several of these means can be embodied by one and the same 
item of hardware. The mere fact that certain measures are recited in mutually different 
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dependent claims does not indicate that a combination of these measures cannot be used to 
advantage. 

In summary, the invention provides coding and decoding of an audio signal 
including estimating a position of a transient signal component in the audio signal, matching 
5 a shape function on the transient signal component in case the transient signal component is 
gradually declining after an initial increase, which shape function has a substantially 
exponential initial behavior and a substantially logarithmic declining behavior; and including 
the position and parameters describing the shape function in an audio stream. 
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1. A method of encoding (1) an audio signal (x), the method comprising the steps 
of: 

estimating (1 10) a position of a transient signal component in the audio signal; 

matching (1 1 1,1 12) a shape function on the transient signal component in case 
5 the transient signal component is gradually declining after an initial increase, which shape 
function has a substantially exponential initial behavior and a substantially logarithmic 
declining behavior; and 

including (1 5) the position and shape parameters describing the shape function 
in an audio stream (AS). 

10 

2. A method as claimed in claim 1 , wherein the shape function is a Laguerre 
function or a generalized discrete Laguerre function. 

3. A method as claimed in claim 2, wherein shape function is a Meixner function 
15 or a Meixner-like function. 

4. A method as claimed in claim 2, wherein at least one of the shape parameters 
is determined by a ratio of slopes of running first and a second order moments of the audio 
signal (x). 

20 

5. A method as claimed in claim 1, wherein the shape parameters include a step 
indication in case the transient signal component is a step-like change in amplitude. 

6. A method as claimed in claim 1, wherein the position of the transient signal 
25 component is a start position. 

7. A method as claimed in claim 1, the method further comprising: 
flattening a part of the audio signal that is furnished to at least one sustained 

coding stage by using the shape function in a gain control mechanism. 
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8. Method of decoding an audio stream, the method comprising the steps of: 
generating (3 1 ) a transient signal component at a given position; and 
calculating (31) a shape function based on received shape parameters, which 

5 shape function has a substantially exponential initial behavior and a substantially logarithmic 
declining behavior. 

9. Audio coder (1), comprising: 

means for estimating (1 10) a position of a transient signal component in the 

10 audio signal; 

means for matching (1 1 1 ,1 12) a shape function on the transient signal 
component in case the transient signal component is gradually declining after an initial 
increase, which shape function has a substantially exponential initial behavior and has a 
substantially logarithmic declining behavior; and 
15 means for including (15) the position and shape parameters describing the 

shape function in an audio stream (AS). 

10. Audio player (3), comprising 

means for generating (31) a transient signal component at a given position; 

20 and 

means for calculating (3 1 ) a shape function based on received shape 
parameters, which shape function has a substantially exponential initial behavior and a 
substantially logarithmic declining behavior. 

25 11. Audio system comprising an audio coder (1) as claimed in claim 9 and an 

audio player (2) as claimed in claim 10. 

12. Audio stream (AS) comprising: 

a position of a transient signal component; and 
30 shape parameters describing an shape function which shape function has a 

substantially exponential initial behavior and a substantially logarithmic declining behavior. 

13. Storage medium (2) on which an audio stream (AS) as claimed in claim 12 has 
been stored. 
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